  x = [0  .901109703749706  .651999996844527  -.639170076615181  .846554698529898  .0345144275346999  -.827655806165909  .291785267894841  .266944316155782  -.284065879033402 -.288072277552419  -.0186527238087105  -.623647682472626 -.915404344049402  .603689131707478  ];

  y = [0  0.  -.758132414106013  .477025636806674  .000406172446199785  -.756421284141003  -.416002940782625  .758519952078692  -.477382835772259  -.941207434336633  .941077403195449  .415901906368019  -.313787114383768  .313689208089463  .756536501280898  ];


  z = [1  -.433591168969247  .0114562994254021  .603248004546166  .532301585025724  .653173480164671  .376732557366718  .582673870821469  -.837165670690192  -.182852787559919 -.177165696738032  -.909218169733735  -.715961741291552 -.25225774046338  -.251419876875175  ];
  
  
    plot3(x,y,z,'rs','Markersize',15);
  sphere;
  axis equal;
  hold on;
  plot3(x,y,z,'rs','Markersize',15,'Markerfacecolor','black');
  hold off;
  
  s = 0;
  
  for i = 1:length(x)
      for j = i+1:length(x)
          s = s + (x(i)-x(j))^2 + (y(i)-y(j))^2 + (z(i)-z(j))^2;
      end
  end
  s